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Abstract 

The bulk viscosity is introduced to model unified dark matter. The viscous unified model assumes 
the universe is filled with a single fluid with the bulk viscosity. We review the general framework 
of the viscous cosmology. The Hubble parameter has a direct connection with the bulk viscosity 
coefficient. For concrete form of the bulk viscosity, the Hubble parameter which has the scaling 
relation with the redshift can be obtained. We discuss two viscosity models and the cosmological 
evolution to which they lead. Using SNe la data, the viscosity model can be fitted. We briefly review 
the fitting method here. 
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I. INTRODUCTION 



Both dark matter problem and the cosmic acceleration problem challenge physicists' un- 
derstanding of the universe. In the standard ACDM model, two mixed fluids, dark matter 
and dark energy fluid, are assumed. These two fluids influence the cosmic evolution sepa- 
rately. However, present gravitational probe does not have the ability to differentiate these 
two fluids. This is the dark degeneracy problem [l| j2|. It is reasonable to model dark matter 
and dark energy with single fluid or single field assumption. Some unified models have been 
proposed to detect the possibility of this unified assumption, like unified dark fluid model 
P ^ jsl sl which assumes the siiigle fluid equation of state; Chaplygin gas model and 

isl - which discuss the cosmology consequences of 



generalized Chaplygin gas jsl 10 1 11| 12 1 



an exotic equate of state;"ca^; Md met^d 000. 

The introduction of viscosity into cosmology has been investigated from different view 
poi nts 17| 181]. There are some recent developments like dark energy model 19|] 20| 21 1 



22|, the cosmic singularity 23|. In this review, we give a brief introduction to unify dark 



matter and dark energy with viscosity medium. In such models, the universe is assumed to 



be filled with viscous single fluid [24 |25| |26| |27| |28| |29( |30|. The cosmic density is not 
separated as dark energy part and dark matter part. The bulk viscosity contributes to the 
cosmic pressure, and plays the role as accelerating the universe. After considering the bulk 
viscosity, the cosmic pressure can be written as 



p=(7-l)p-3Ci/ 



(1) 



Where 7 parameterizes the equation of state. Generally the form of bulk viscosity is chosen 
as a time-dependent function. In [3l| [32| [33| [3J], a density-dependent viscosity ( = ap^ 
coefficient is investigated extensively. For modeling the unified dark matter and dark energy, 
it is often assumed that the parameter 7 = 1, that the pressure of the viscosity fluid is 
zero and the viscosity term contributes an effective pressure. There raises some problems 
here. From the observational results 35|, the cosmic density nearly equals to the cosmic 
pressure. In the viscosity model, the viscosity term dominates the cosmic pressure, and 
surpasses the pressure contributions from other cosmic matter constitutions, which contradicts 
the traditional fluid theory. 36|] 37| propose non-standard interaction mechanism to solve 
this problem. Obviously, it is important to build solid foundation for the research of the 
viscous cosmology. 

Equation of state w < —1 lies in the phantom region. It is shown that cosmology mode 
with such equation of state possess the so-called the future singularity called the Big Rip 



The larger viscosity model parameter space can help to solve the cosmic singularity problem 
and produces different kinds of evolution mode of the future universe, for more details 29|. 

The rest of this review is organized as follows: In the next chapter, general framework 
of the viscosity model will be reviewed. In Sec. Ill, we discuss the modeling of the unified 
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model with viscosity. In this section, two concrete models are analyzed. In Sec. IV, data 
fitting method is introduced briefly. 



II. GENERAL FRAMEWORK 

We consider the standard Priedmann-Robertson- Walker metric, 

ds' ^ _dt2 + a{tf{-^^^ + r'dn') . (2) 

For the sake of simplicity, we choose the flat geometry k — 0, which is also favored by the 
update result of the cosmic background radiation measurement. 
The general stress-energy-momentum tensor is 

= (P + P)U^U, + pg^, - V> (3) 

where ( is the bulk viscosity. The expansion factor 9 is defined by 9 — U.^^ — 3^, and 
the projection tensor /i^,^ = Qi^^, + U^U,^. In the co-moving coordinates, the four velocity 
C/^ = (1, 0, 0, 0). We do not specify the concrete form of ( in this section. Generally speaking, 
C is a quantity evolving with time t or the scale factor a{t). We will see below that non-trivial 
and more complicated ( can produce different results especially useful for the late universe 
modeling. 

From the usual Einstein equation, 

C^. = i?/.. - ^Rgf^u = SnGT^,, (4) 



we obtain two equations which we call the modified Friedmann equations: 



(5a) 



3 

d 4:TtG , . , , 

- = ^P + 3p, 5b 

a 3 

where p is an effective pressure, p = p — C,9. 

The covariant conservation equation T^^ = yields 

p+(p + p)^ = 0. (6) 

The existence of a bulk viscosity contributes a modification to the pressure p, thus we see 
the Friedmann equation and the covariant conservation equation are invariant under the 
transformation 

P^P^P-CO, (7) 
The covariant energy conservation equation becomes 

p+{p- (9)9^0. (8) 
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If define the dimensionless Hubble parameter liere 



H^_ p 

t2 



Pc 



(9) 



3H 

where pcr = g:;;^ is the critical density now. Using the dimensionless Hubble parameter, Eq. 
dH]) can be transformed as 

1 d(/l2^ 



+ 3h-' = 9Xh^ 



(10) 



Ho dt 

where the bulk viscosity is redefined as A = y^. Through the simple relation between scale 
factor a{t) and the redshift z 

dt = -^da, (11) 
aH 

we transform Eq. fIlOp into a differential equation with respect to the scale factor a{t) 

dH 3 3C 

— H H = — . 

da 2a 2a 

Solving this equation, we obtain a integral form of H{a) 



(12) 



H(a) =Cia~^/^+ / — exp( / —da) da exp( 
'2a ^ / / 



2a 



2a 



da) . 



(13) 



Different forms of viscosity can be used here to make this integral calculable, numerically or 
exactly. 



III. UNIFIED SINGLE FLUID 
A. Redshift-dependent model 



In [39[, authors assume the bulk viscosity takes the form as an Hubble parameter depen- 
dent function. A redshift-dependent viscosity is proposed in [3]. This bulk viscosity is a 
combination of a constant and a scaling relation term 

9A = Ao + Ai(l + 2)", (14) 

where n is an integer, Aq and Ai are two constants, which will be fitted from the observational 
data sets. 

After taking account of this ansatz, the integration is easily to work out. We get 

2V ^ ; ^3 2v ^ ; 3(2n-3)^ ^ (2n-3)2^ ^ 2n - 3^ ^ 9 

(15) 



4 



Since we have assumed the spatial flat of the universe, the consistency condition requires 
h{0) = 1. Thus this sets a constraint on the model parameters as 

^ = 1-A, + ^^ (16) 
3 ^ 2n-3 ^ ' 

We remind the readers that we make the single fluid assumption above, and we do not 

concretely specify the constitutions of the cosmic density p. In this single fluid model, values 

of model parameters Aq, Ai and A2 will be fitted, and their meaning are not explained. But 

when we compare it with two-fluid model, that the universe is filled with dark matter and 

dark energy fluid, more constraints can be added. The solution consists terms with different 

scaling relation. The first term has the form like C(l + z)^, which have the same evolution 

behavior as the cold dark matter. Their simplicity leads us to correspond parameter A2 to 

dark matter ratio Qm 

Xl = ^m. (17) 

This identity can help us utilize more data to constrain the viscosity model. The result is 
also consistent with that obtained from the standard model(ACDM). The shift parameter TZ 



40( | 4l| and the distance parameter A is defined as 



§-y (18) 



and 

= \/ 1^772 



respectively. Both of them are dependent on dark matter ratio Qm, and in the joint statistical 
analysis they provide strong constraint on Q^n- 



B. Effective equation of state model 



Another viscosity model reviewed here is proposed in 26|, where a general form time- 
dependent viscosity is discussed 

C = Co + Ci- + C2-. (20) 

a a 

An interesting feature of this model is its effective equivalence to the following equation of 
state 

p = (7 - l)p + Po + whH + WH2H^ + WduH (21) 

where po, wh2 and Wdh are free parameters. The corresponding between two groups of 
coefficients are 

wh = -3Co, (22a) 
WH, = -3(Ci + C2), (22b) 
WdH = -3C2. (22c) 
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The parameterized bulk viscosity combines terms related to the "velocity" a and "accelera- 
tion" a, which can be seen to describe the dynamics of the cosmic non-perfect fluid. After 
eliminating p and p, a differential equation about the scale factor a{t) can be obtained 



1 {K-^)WdH 



/U/x Z 



+ 



+ {Ky2)WdHa l+{Ky2)WdH 



(23) 



Another feature of this model is that this differential equation can be solved exactly, and the 
evolution function of the scale factor a{t) is definite. This evolution function is especially 
convenient for discussing the cosmic singularity. 

With the initial conditions a{to) = Oq and 9{to) = Oq, when 7 7^ 0, the scale factor can be 
obtained as 



ait) = «o{^(l + 7^or-^)exp[^(i + l)] + 



i(l-7^oT+|)exp[-^(l 



T 
1 

f 



(24) 



where the parameters are redefined as 
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7 + {Ky3)wH2 

i + {Ky2)waH 

' l + (Ky2)WdH 

-(«V2)P0 



1 

1 _ 

1 



T? 



1 67 
7? + If- 



Prom Priedmann equation, p can be written as 



Pit)- 



3k^7^ 



:i + 7^oT-f)(^ + ^)exp( 



t-to ■ 
T ' 



:i-7^oT+f)(i-^) 



:i+7^oT-^)exp(^^: 



(l-7^^oT+^) 



(25) 

(26) 
(27) 
(28) 

(29) 



The model parameters leave enough space to produce various evolution behavior, which can 
be interpreted in different ways. In this review, we emphasis its power to unify dark energy 
and dark matter with the single fluid assumption. According to the parameters redefined 
above and the Friedmann equation, the equation of state can be converted to 



p = (7 - l)p 



(30) 



The case 7 = and Ti — > 00 corresponds to the ACDM. With the aim to unify dark energy 
and dark matter, the case 7 = 1 and 72 00 is especially considered. This case corresponds 
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to a single fluid with constant viscosity. The relation between p and p can be obtained from 
the general equation of state above 

Therefore, it is straightforward to eliminate p from the covariant energy conservation equation, 
and to work out the solution of p. Using Friedmann equation, H{z) can be obtained 

H{z) = Ho[n.Xl + zf^ + {l-n.,)]. (32) 

VL^ is the only one model parameter. Its value can be fitted from SNe la observational data. 



IV. DATA FITTING 



We review the method to fit the model parameters. More details are illustrated in 42 1. 
The data sets we use are SNe la, BAO and CMB. The 397 Constitution sample j43| combines 
the Union sample 44 1 and the low redshift {z < 0.08) sample j45|. The co-moving distance 
cIm in FRW coordinate is 



d 



M 



Hiz' 



-dz' 



The apparent magnitude which is measured is 

m = M + 5\og^oDLiz) 
where the dimensionless luminosity = H^di^z) and 

dL = {1 + z)dM{z). 



(33) 



(34) 



(35) 



where M is the absolute magnitude which is believed to be constant for all SNe la. In the 
SNe la samples, data are given in terms of the distance modulus pobs = rni^z) — Mobs{z). The 
for this procedure is written as 

> n 2 



X 



E 

i=l 



f^obs{Zi) — Pth{Zi] Co) 



^obs \Zi 



(36) 



where pth means the distance modulus calculated from model with parameters Cq (a = 
0,1,2...). Together with the shift parameter TZ and the distance A, the total xlotai 
joint data analysis is 



Xtotal 



n-n 



obs 



A- A 



obs 



(37) 



For the redshift-dependent model, the relation between distance modulus and redshift is 
plotted in FIG.l. The model calculated value and the Constitution data is compared in the 
figure. 
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FIG. 1: Relation between distance modulus and redshift. The solid line corresponds to the theo- 
retical value calculated from model concerned. The dots with error bar are the data from the 397 
Constitution sample. 




V. CONCLUSION 



In this review, we discuss three aspects of the viscosity model, 



General framework for viscosity modeling. General form of Hubble parameter is pre- 
sented. This general form is convenient for comparing different scale factor(or redshift) 
dependent viscosity models. 

Two kinds of viscosity models are used to model unified models. 

Observation constraint is necessary for model building. We can see the fitting results 
are consistent with data. It is prospected that more accurate direct measurements of 
Hubble constant will provide a new constraint on cosmological parameters 46 1. 



Especially we focus on its application on modeling the unified dark energy and dark matter. 

In the cosmic background level, dynamical analysis can be performed. The statefinder 
method is useful for discriminating different models 47(] 48|] 49|] 50|] 5l|- Compared with 
ACDM model, evolution of the statefinder of the viscosity model is different and can be 
discriminated easily, more details can be found in 52| 29|]. More plentiful and accurate data 
will improve the power of the statefinder method, which will give enough constraint on the 
late universe model. 

We review the viscosity model which is on the level of zero order. The perturbation 
analysis and the large scale structure are especially useful for the model building. The model 
predictions need to be consistent with CMB and LSS data. Some works has investigated 
the perturbation aspects of the viscosity model 53|] 3J]. After corresponding the model 
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parameters, the viscosity model has the connection with the Chaplygin gas modeL Though 
the Chaplygin gas model can fit the SNe la data well, in the perturbation level it is found the 
Chaplygin gas model does not behave in a satisfactory way. Whether the viscosity models 
could behave well needs further investigation. 
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